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H- ■ 

^y. , We derive generalized quasiclassical transport equations which include the terms responsible for 

Qs ' the Hall Effect in the vortex state of a clean type-II superconductor, and calculate the conductivity 

tensor for an s-wave superconductor in the high-field regime. We find that below the superconducting 
transition the contribution to the transverse conductivity due to dynamical fiuctuations of the order 
flj ' parameter is compensated by the modification of the quasiparticle contribution. In this regime 

^\ , the nonlinear behaviour of the Hall angle is governed by the change in the effective quasiparticle 

^^ ■ scattering rate due to the reduction in the density of states at the Fermi level. The connection with 

CO ' experimental results is discussed. 

(N 

o 
o 

I ■ I. INTRODUCTION 

5-1 ' 

Oh; 

^ I In the recent years a significant body of work has been devoted to the better understanding of the Hall ejfect in 

the mixed state Cj£ptype-H superconductors, which has remained a theoretical puzzle for almost thirty ycarsy'u The 

phenjamenologicaluij theories predict that the Hall angle in the flux-flow regime is either identical to that in the normal 

stateH or constanta, and the underlying microscopic basis for recent generalizationaa is not well understood. Theories 

which make use of the time-dependeBt Ginzburg-Landau equations (TDGL) find that the Hall conductivity is not 

'"^ modified in the superconducting stateEI. These predictions are at variance with the strongly nonlinear behaviour (3*r* 

^ function of magnetic field) found in experiments performed on both low- Tc materialsB^El and the high-Tc cupratesBfl. 

O For dirty superconductors (Z <C Co- where I is the mean free path and ^o is the superconducting coherence length), 

1^1 - transport coefficients can be determined from microscopic theory by a straightfocjaiard expansion in powers of the 

order parameter, A. The results of such a calculation for the transverse resistivityQ'u explain qualitatively the sharp 

T-H , increase in the Hall angle below the transition observed in experiment (although, to our knowledge, no systematic 

^ ' comparison has been made), and provide the physical basis for a generalized TDGL approach, in which the relaxation 

^^ , rate is assumal. to be complex, rather than purely real, to allow for a modification of the transverse transport 

^-j ■ coefficientsld'LJ'Eil. The small parameter in the expansion of the microscopic equations is proportional to both the 

^vj order parameter and the mean free path, therefore, it is not small in the cleaii.II ^ ^q) limit. In this regime a 

y—( straightforward expansion is not possible; the TDGL equations are not applicablell3li3, and so an alternative approach 

f — . ' is needed to determine the transverse transport coefficients. 

0^ ■ In this work we develop an approach to calculate the transport coefficients, including the Hall Effect, of clean 

.<— > ■ type-II superconductors in the vortex state and present the results of a calculation of the Hall conductivity of a clean 

2 r s-wave superconductor in the mixed state near the upper critical field, i?c2- The method is based on the quasiclassical 

approximation to the microscopic theory, due originally, in the context of superconductivity, to EilenbergerEfl and 

Larkin and Ovchinnikovlla, which we generalize to include the terms responsible for the Hall Effect in a charged 

superfluid. We solve the equations of this quasiclassical theory to obtain the longitudinal and transverse resistivities 

Q in the mixed state. We choose to consider an s-wave superconductor as both the normal state and superconducting 

CJ . properties of the low- Tc compounds are well known, and comparison between theory and experiment is fraught with 

sj ' less ambiguity; however the approach developed here can easily be generalized to consider superconductors with other 

• i-H . than s-wave symmetry. 

p\^ • The microscopic Green's function contains all the information about the single particle properties of the system. 

JH [ In particular, it oscillates on length scales of order of the inverse Fermi wave vector kj . However, when calculating 
. - - 1 transport coefficients, we are for the most part only interested in the long-wavelength response. It is then sufficient 
to determine the envelope of the Green's function rather than its detailed form. In the quasiclassical approach the 
rapid oscillations associated with the presence of the Fermi surface are integrated out of the basic equations and 
slower varying quantities such as external fields or the self energy are expanded around their values at the Fermi 
surface. The resulting transport like equations contain the microscopic physics relevant to the problem and are easier 
to solve. The basic premise of quasiclassical transport theory is that all macroscopic physical quantities vary slowly on 
a microscopic length scale, and that all the relevant momenta are small compared to the Fpmi momentum pf . This, 
approximation has been applied successfully to study transport phenomena in superfluidstSI and superconductorsEZI 



and to investigate the behavior of the unconventional superconductorgla. Recently it has been usedpto analyse the 
most relevant contributions to the Hall effect in a dirty supercondustor in the limit of isolated vorticesll3 as well as to 
investigate the forces acting on a single vortex in the clean regimaH3. 

In the next two sections we present a derivation of the generalized quasiclassical equations, which include all 
the terms contributing to the Hall Effect in the mixed state of a clean type-H superconductor in the high-field 
regime. Section || introduces a general quasiclassical formalism and the basic ideas involved in the analysis of 
transverse transport in the quasiclassical approximation, illustrated by application to the simple case of a normal 
metal. We show how the standard Drude results for longitudinal and transverse conductivity are obtained within 



this quasiclassical approximation. In Section HI we use the same approach to derive a generalization of the standard 
quasiclassical approximation for superconductors to include the terms responsible for the transverse conductivity and 
obtain linearized transport like equations for a clean superconductor. To soUk this system of equations near the 
upper critical field we employ the approximation of Brandt, Pesch and TewordtcJ (BPT), in which the normal part of 
the matrix propagator is replaced by its spatial average over a unit cell of the vortex lattice, while the exact spatial 
dependence of the order parameter is retained. Using an operator formalism, we are able to solve the leading order 



equations for the distribution function in Section [V , and obtain the longitudinal and transverse conductivities within 



linear response theory in Section M and Section VI respectively. In the last section we summarize the results and 
compare them with the existing experimental data. 

II. QUASICLASSICAL APPROACH TO TRANSPORT IN A NORMAL METAL 
A. Mixed Representation and the Standard Quasiclassical Equations 

Our starting point is the microscopic Dyson's equation 



-|:-CHVx)-yAs(x, 



y) 



G{y,x') = 5{x-x') (1) 



for the Green's function 



G{x,x') = -{Tri^{x)^\x')). (2) 



Here '4'{x) ^nd ^^^(2;) are field creation and annihilation operators, which depend on the four vector x — (x, r), 
angular brackets denote the statistical average, and the operator Tr arranges the field operators in ascending order 
of imaginary time r. In Eqn. ^ C is the single particle energy operator, and S is the self energy which may be due 
to interactions or impurity scattering, its exact form has to be determined from microscopic considerations. Dyson's 
equation can also be written in the form 



G{x,x') 



d 



d 

7 - C(+*Vx' 



- d'^yd{x,y)Y.{y,x') = S{x-x'), (3) 



The operators in this equation are understood to act on the Green's function on their left. It should be emphasized 
that Eqs. (m and (ra) contain the same physical information and only differ in the form of writing, i.e. the same 
function G satisfies both. We will use the terms right-hand and left-hand Dyson's equation for Eqs. (|l|) and (|^) 
respectively. l_. _ 

The derivation ot-Lhe quasiclassical equations given here follows the general approach of Rainer and SerenetSS and 
Eckern and SchmidEj. First we consider the linear response of a metal to a constant uniform electric field described by 
a vector potential A(t) = Aexp(ia;oT). To incorporate the vector potential into the microscopic equations we replace 
the momentum operator by its gauge invariant counterpart C(~*Vx) ^ C(^*^x — ^-^.(t)), and expand this expression 
to obtain terms linear in the external field. To integrate out the rapid oscillations associated with the presence of the 
Fermi surface we first change variables from x and x' to center of mass and relative coordinates R = (x + x')/2 and 
r = X — x', and carry out a Fourier transformation in the latter according to 

G(p, R)^J d'rG{R + |, R - ^) exp(-zpr) (4) 

In a translationally invariant system the Green's function only depends on the relative coordinate. Therefore, depen- 
dence on the position of the center of mass R appears only in the presence of external fields. To treat the effect of 
slowly varying fields quasiclassically we expand in quantities varying on the length scale of the wavelength of these 



fields, which is equivalent to expanding in powers of Vr. If A{k, — iVx) is a local operator which depends only on 
position and momentum and acts on the Green's function G{x,x'), then 



d^r exp(-ipr)v4(x, -iVx)G(a;, x') = / dhc^p{-ipr)A{R + |, -iVr - ^Vr)G(R +^^'^-1) 
= y d-\A(R+ ^Vp,p- ^Vr)G(R+ |,R- ^) expHpr) = A(R+ ^Vp,p- ^Vr)G(p,R). 
The final expression can be written as A o G, where the "circle-product" is defined astJO 



(5) 



^(p,R)oB(p,R) =exp(-(Vp,VR, - Vp,VrJ)v4(pi,Ri)B(p2,R2)|ri=r,=r, 



(6) 



Using this definition the right and left-hand Dyson's equations can be written in terms of the mixed set of variables 
p and R as 



d 



dr 
G(p,R;r,r') 



C(p-eA(r)) 
d 



dr' 



oG(p,R;T,T')-y"dTiS(p,R;T,ri)oG(p,R;Ti,r')-<5(T-r') 
C(p - eA(T')) - y dTiG(p,R;T,Ti) o I](p,R; ri,T') = 6{t-t'). 



(7) 
(8) 



Direct expansion of Eqs. (0) and (g) in powers of the spatial gradient is not possible since in the definition Eq. (g) of 
the circle-product this gradient is coupled to derivatives with respect to momentum, and the Green's function varies 
rapidly with momentum near pf. To avoid this problem we make a transformation from the set of variables (p, R) to 
the set (s, C, R), where s is a parameterization of the Fermi surface, and integrate the quantum mechanical equations 
over the quasiparticle energy C before expanding. The integrated Green's function 



g(s,R;r,r') = ^ydCG(p,R;T,T') 



(9) 



only depends on the components of momentum parallel to the Fermi surface and the remaining dependence on p and 
R is slow. We now transform Eqs. (0) and (||) for the full microscopic Green's function G into equations for the 
quasiclassical propagator g. This quasiclassical propagator will play the role of a distribution function in the resulting 
transport like equation. 

Let us first compare terms of zeroth order in the gradient expansion of Eqs. (0)and (^. Since the imaginary 
time T varies between and l/T", where T is temperature, the first term in the equation gives, after integration, a 
contribution of order Tg. If we assume that the self energy varies slowly for momenta close to the Fermi momentum 
IpI ~ Pf, we can approximate 



y dCS(p, . . .)G(p, . . .) « I](p/, . . .) y dCG(p, . . .) = a{i 



.)g{s,---). 



(10) 



On the other hand, the term involving C gives a much larger contribution since the integration region includes C ~ e/. 
Because of this term and the delta function on the right-hand side the equations cannot be integrated directly. Instead, 
we subtract Eq. dsh from Eq. (R) to obtain a homogeneous form before integrating term by term and expanding in 
the gradients. The zeroth order term involving CG then cancels. Expanding to first order we obtain 



y dC[G o C(p - eA(T')) - C(P - eA(T)) o G] 
where the Fermi velocity is defined as 

dC 



A(r')-A(r) 



g + jvVr5, 



dp 



(P/)- 



(11) 



(12) 



If the spatial dependence of the distribution function is determined by the wave vector q of an external field, the 
product vq is not necessarily small compared to the temperature and the self energy, so that this term has to be 
retained in the leading order equation. Since the small parameters in the expansion are of order l/[kfX), where A is 
a typical wavelength of the electric field, for the terms involving the external vector potential, or, if the self energy 
is due to impurity scattering, l/{kfl), C must always be expanded to one order higher in small quantities than other 



terms in order to obtain a contribution of similar order. It should also be emphasized that, since there are several 
small parameters in the problem, it may be necessary to expand terms to different order in gradients to account for 
all the contributions to a particular physical quantity. 

It is convenient to Fourier decompose the integrated Green's function into Matsubara frequencies 



g{s, R; r, t') ^T^ g{s, R; a;„, a;„') exp(-iw„T + iw„'r'). 



(13) 



where cOn are the fermionic frequencies LUn = (2n + 1)ttT. Then the basic transport equation of the quasiclassical 
formalism becomes 



iuJn - iuJn' + «v(s)Vr 



9[s 



R;a;„,a;„/) + evA 



g(s, R; ujn - uq, w„') - g{s, R; a;„, a;„/ + loq) 



-T 



E 



t(s, R; cj„, U!k)gis, R; Wfc, w„') - g{s, R; uJn, ujk)a{s, R; ujk, w„') 



= 0. 



(14) 



The exact form of the self energy, a, is determined from microscopic theory. In principle, all higher order terms in 
the spatial gradient can be included in this equation consistently using the definition of the circle product. 

It should be noted that, in the absence of a perturbing potential, or impurity scattering leading to the appearance 
of the self energy, the Green's function is independent of the coordinate R and is diagonal in frequency space, and 
therefore Eq. ( |l4| ) is trivially satisfied by any function g. This is not surprising since in subtracting the right 
hand Dyson's equation from the left hand equation the information about a particular solution of the inhomogeneous 
equation has been lost. The particular solution describes the unperturbed non- interacting electron gas, and is obtained 
by integrating the function Go = [iuJn — C(p)]~^ over the quasiparticle energy to find the quasiclassical distribution 
function of a normal metal, go = —isgn{uJn) ■ This function serves as input for any perturbative approach to transport 
in a metal. 



B. Semiclassical Treatment of the Magnetic Field and the Lorentz Force. 



Eq. (nj) is sufficient to analyse longitudinal transport in a normal metal but it has to be generalized to determine 
the Hall conductivity. If the vector potential ^(R) describing the magnetic field is taken to be of order "smalF , the 
field itself, H = V x ^, becomes of order "(small)'^" and the Lorentz force, which is proportional to both electric 
and magnetic fields, disappears from the perturbative expansion of the quasiclassical equations. This observation led 
Rainercj to point out that in order to analyse the Hall Effect in a normal metal, the vector potential ACR) must 
be considered as a leading order quantity and should be included in the equations semiclassically rather than being 
treated perturbatively. Now the quasiparticle energy ( depends on the generalized momentum p — e A — ey^(R) . This 
replacement is exact. The semiclassical approximation, which is applicable in the long wavelength limit where the 
quasiclassical approach is appropriate, treats the momentum operator as a c-number. Therefore in the transformations 
described in Eq. (|5|) the momentum p and the coordinate R are no longer independent variables, rather, they are 
coupled by the presence of the vector potential, which depends upon the coordinates. As a result the gradient expansion 
of the integrated Green's function cannot be carried out independently in the Fermi surface parameterization s and 
the spatial variable R. For a general transformation of variables from the set (p, R) to the set (C, Si, R) 



d 



9C 9 dsi d 



d 

dRa dRa dRa d( ' dRa dSi 

d dQ d dsi d 

dpa dpa dC dpa Osi 



(15) 
(16) 



where the derivatives on the right hand side are computed at constant CjS,R rather than p,R. Using the explicit 
semiclassical R-dependence of C and s 



= — e- 



dR. 
ds, 



dpp dRa 

_ dsi dAf3 

dRa dpp dRa ' 

we obtain from the expansion of the terms involving ( 



-ev/B 



dA, 



dRa' 



(17) 
(18) 



dC 



-C(p - eA(T) - eAiR)) o G + G o C(p - eA(r') - e^(R)) 



(19) 



/ dC ( -C(P ~ eA(T) - eA(R))G + C(p - eA(r') - eyl(R))G 



+ i 



dC dG dC dG 
9p 9R 9R 9p 



A(r')-A(r) 



9 + «vVR.g + ie(v x H)- — , 
CP|| 



where p|| denotes the component of the momentum p parallel to the Fermi surface. In the last line of Eq.(n9h we have 
used the result 



dC dG dC dG 
dpdn^ffR.'dp 



dC dG 



dpa dRa 

dG 



dC dsi d( dsi 



dRa, 



+ eva 



dpa dRa dRa dp, 
dAa dAp 



dG 

dsi 



dRp dRa 



ds, dG „ ^ , „, dG 

TT^T^ = vVrG + e V X n)^—. 

dpp dsi dp[| 



(20) 



The new term is the familiar Lorentz force driving term of the classical Boltzmann transport equation. Here it appears 
from taking into account correctly the semiclassical dependence of the momentum on the external field. The basic 
quasiclassical equation Eqn. (|lj) now takes the form 

g(s,R;u;„,w„') + 



iv(s)VR + ze(v X H) 



^Pll 



lUJn — li^n 



5(s,R;w„,w„ 



+evA 



cr(s, R; w„, ujk)gis, R; ujk,ujn') - g(s, R; w„, Wfc)cr(s, R; ujk, ujn') 



(21) 



C. Linear Response 

In general Eq. ( pl| ) is a nonlinear equation. To calculate transport coefficients it is sufficient to keep only the terms 
linear in the external perturbation - in this case in the electric field - and determine the Green's function g within 
linear response. We separate the propagator into a leading term and a part linear in the perturbing potential 



g = goii^n)7f;S^„,u„ 



.gf^) (s, R; w„, cjo) -<5^„,^„,+u;o ■ 



(22) 



If the self energy is due to elastic impurity scattering, it can be separated in a similar way into ctq Eind a^^' . As noted 
the equation for the leading order terms go and ao is satisfied trivially; the terms of linear order are given by 



iujQ + o-o(-) - ctq + iv(s)VR 



,(1) 



+ ie(v X H) 



dg^ 

dp\\ 



(evA + aW)(5o( 



go 



(23) 



here we have used a shorthand notation go = 50 (^n) and go(~) = goi^n — ^o)- This equation is the basis for the 
analysis of transport in a normal metal. It has to be solved together with the self consistency condition relating the 
change in the self energy to the modification of the Green's function g^^\ 

Since, throughout this work, we will be concerned_with the electrical conductivity, we have to define the current in 
terms of the distribution function. It is well knowntj that, if in the microscopic equation for the current density 



j(x) = -r5] /d3ppGW(p,R^x;. 
m ^-^ J 



Ne' 



(24) 



the integration over energy is carried out before summing over frequencies, the contribution from the high energy 
regions (far above and below the Fermi surfaj^ej) exactly cancels the diamagnetic term in Eq. (|4|). Then the quasi- 
classical expression for the current becomesc 



i{R)^nN{0)eTj2 f d^s^r{s)g'■^\s,R; 



UJ 



(25) 



where A^(0) is the density of states at the Fermi surface. The problem of calculating the transport coefficients of a 
normal metal is now fully defined. 



D. Conductivity of a Normal Metal 

As an example of the usefulness of the quasiclassical method we will use it to determine the conductivity tensor of 
a normal metal in a magnetic field. We consider an experimental arrangement with constant electric and magnetic 
fields E = Ex and H = Hz. We also assume a spherical Fermi surface 

V = t;(sin6'cos0, sin0sin0, cos^) (26) 

and include the effect of isotropic impurity scattering in the Born approximation, so that the self energy is given by 

a^^jcPsg (27) 

where r is the quasiparticle lifetime. The unperturbed Green's function is given by g^ — — isgn(w„), and therefore 
(To = -isgn(a;„)/2T. 

First consider the longitudinal dc conductivity. In the absence of a magnetic field Eq. (p3), becomes 

iujo + — (sgn(a;„) - sgn(w„ - ujq)) g^^'> = -i(evA + cr^^^) (sgn(a;„) - sgn(w„ - ujq)). (28) 

Since the driving term in Eq. (Eq) is proportional to vA, it is evident that the angular dependence of g^^' is given 
by that dot product, and there is no correction to the self energy since the angular average of 5'-^-' vanishes. Then 
it is obvious from Eq. ( p8|) that g^^' = when w„ and a;„ — wq have the same sign. Otherwise, in the intermediate 
frequency region where Wo > a;„ > 0, 

5« = -'-^. (29) 

UJQ + l/T 

Integrating over the Fermi surface, carrying out the summation in the definition of current density, and analytically 
continuing to the real external frequency according to iluq — > o) + i(5, in the dc-limit (cD — > 0) we recover from this 
solution the standard Drude theory result for the current 

j = ^N{0)eWE = a„E. (30) 

We now turn on the magnetic field. Writing the expression for the Lorentz force in spherical coordinates it is easy 
to check that 

and the linearized transport equation becomes 

i d 

[iujo + — (sgn(a;„) - sgn(w„ - ujo)) - iujc-^]g'^^^ = -iev(s)A(sgn(a;„) - sgn(w„ - ojq)). (32) 

Again, the response function is non-zero in the intermediate region only. In the regime uJcT ^ 1 it is sufficient to solve 
the equation perturbatively, namely, 

(1) 2ev(s)A 

9h = — ——riz + Sg (33) 

UjQ + l/T 

2e2 

Sg = r— :-TruJcV A sin 9 sin (j). (34) 

[ujQ + I/tY 

The transverse current obtained from the correction 5g is, as expected, 

jy = -OniOcTE. (35) 

We have therefore reproduced the results of the Drude theory using the quasiclassical formalism. 



III. QUASICLASSICAL EQUATIONS FOR A SUPERCONDUCTOR 

In this section we generalize the approach developed in Section |l^ to derive a set of quasiclassical equations which 
can be used to analyse both longitudinal and transverse transport in superconductors. 

A. Gorkov equations 

Gorkov's equationsEH for a matrix Green's function G replace Dyson's equations in a fully microscopic approach to 
a superconductor. The diagonal elements of the matrix Green's function 



G = 



G -F 

F^ G 



are the particle and hole propagators, 



G{x,x') = G{x\x)] 



(36) 



(37) 



for singlet pairing, the off diagonal elements are related to the probability amplitudes for the destruction or creation 
of a Cooper pair by 



{idy)al3F(x,x') = -{Tr'lj)a{x)'lj)l3{x')) 

iidy)^pF^{x,x') = {T,i:i{x)^l{x')) 
where ay is the Pauli matrix. Then the right and left-hand Gorkov equations are 



d - ' 

^dz - C(-iVxCTz) + A(a;) 



G(a;, x) — Id yT,{x, y)G{y, x') = 5{x — x')! 



G{x,x') 
The matrix order parameter 



-TTl^z - C(+*Vx'az) + A(x') 



- d^yG{x,y)i:{y,x')^S{x^x')l. 



A = 



A 



-A* 

is related to the off-diagonal elements of the Green's function by 

A{x) =gF{x + 0,x) 
A*{x) =gF\x + 0,x), 

where g is the coupling constant. 



(38) 
(39) 



(40) 
(41) 

(42) 



(43) 
(44) 



B. Quasiclassical Approximation 



The general approach to the derivation of the quasiclassical equations for superconductors is exactly the same as 
that of Section O. We introduce the vector potentials A(t) and A{x) of an electric and magnetic field into the energy 
operator, transform the equations to a set of "mixed" variables p and R by performing a Fourier transform in the 
relative coordinate, and expand in gradients with respect to the center of mass coordinate, after integration over the 
quasiparticle energy. 

Expanding the circle product C(~*Vx) o G to first order in gradients, we obtain 



d^rexp(-ipr)C(-iVx - A{x.) - A)G 



C(p)-2v(VR-2ze^(R) 



vA- — (v X H)-— 
2 ^ ^ dp 



G(p,R). (45) 



On the other hand on expanding the operator C(+*Vx') the combination p -t- zVr. rather than p — jVr, appears after 
Fourier transform in r. Consequently the magnetic field dependent terms arising from the expansion of Vr in Eq. 
( p^ ) have the opposite sign, and 



d3rexp(-ipr)C(+iVx' - yi(x') - A)G 



C(P) + t^v(Vr + 2ieAiR)) - vA(t') 



— V X H h^ 



G(p,R). (46) 



Subtracting Eq. (^6|) from equation (^5|) we regain the result of the Section pB| . The vector potential A appears in 
the expansions in different gauge invariant combinations. This can be easily understood if we remember that operator 
C(— iVx) acts on the annihilation operator ip while operator C(+*Vx') acts on the creation operator ip'' . Then the 
time evolution of the operators describes the motion of particles and holes respectively, and the appropriate gauge 
invariant derivative is different in each case. To determine the transverse conductivity all contributions of the order of 
the cyclotron frequency lOc — eH/mc have to be included in the equations. In a type-II superconductor in the vortex 
state the coherence length ^o sets the length scale for spatial change of the order parameter. Near the upper critical 
field Hc2 the magnetic length A = {2eH)~^/^ ~ ^Oi this immediately implies that the expansion of the operator ^(p) 
has to be carried out not to first, but to second order in spatial derivatives. The second order derivative of C with 
respect to momentum is, by definition, the inverse effective mass tensor, which in the simple case of a spherical Fermi 
surface becomes equal to the inverse effective mass m. In the expansion this term is coupled to square of the spatial 
gradient, so that its contribution 



dpadpi3 



"^c^pG 



1 



lA2 



G oc ujcG 



(47) 



is comparable to that of the Lorentz force term and has to be taken into account. Neglecting terms quadratic in the 
electric and magnetic fields and assuming a Fermi surface with the reflection symmetry C(p) = C(~P) we obtain the 
expansion of the quasiparticle energy operator 



C(-»Vx) - C(P) - ^v(V - 2ieA{K)) vA(r) _ | [v x H] ^ 



-l-{V-2^eAmf 



le 
2m 



A(V - 2ieA{'R.)) 



C(+*Vx) -^ C(P) 



^v(V + 2.e^(R)) + vA(t) + ^ [v x H] ^ 



(48) 



(49) 



^^N + 2ieA{Ii)f 
am ' 



2m 



A(V + 2ieA(Il)) 



and similar expressions for the operators C(+*Vx') and C("*Vx')- 

Now consider the remaining terms in the expansion of the microscopic Eqs. (Hfl) and (ETI). Here we are concerned 
with the change in the Hall conductivity of a superconductor relative to the normal state value. This change involves 
the magnitude of the superconducting order parameter A, which appears in our analysis in the dimensionless combi- 
nation (AA/u). It is then easily seen that linear terms in the gradient expansion of the order parameter have to be 
retained in the equation since a typical term in the expansion 



dAdG A 1 ;:, 

— — -— — oc G : 

oR op A m,v 
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will contribute significantly to the change of transverse conductivity upon entering the superconducting state, 
panding to first order in the gradients we obtain from Eq. (M) 



(50) 



Ex- 



A{x)G{x,x') - G{x,x')A{x') 



A(R, t)G(p, R) - G(p, R)A(R, r') 
ildAdG dddA] i [ dA dG 



2 [aR c»p dp dR.\ 2 [ ap 9R ^R ^p 



dGdA 



(51) 



and, similarly, 



A -^(x, y)Giy, x') + G(x, y^y, x') 



/"drWS(p,R;r,ri)G(p,R;ri,r')-G(p,R;T,Ti)E(p,R,Ti,r') 



9S dG dG as ■ 


i 


'dY.dG dGdY. 


dR dp dp dR 


2 


dpdR dRdp 



(52) 



Using the results of Eqs. (|4§|)-([49|), ( |5l| ) and (|52|), subtracting the left-hand Gorkov-Dyson equation from the 
right-hand equation, and integrating over the quasiparticle energy, we obtain the quasiclassical transport equation for 
a superconductor, which can be written using the matrix notation as follows 

iuJnazd-iuJn'g^z + Ag-gA + ivVg + evAlazd-g^z] - T^—Aldz'^d+^d^z] 
zc (J r 

-I-— (v X H) [dzd + gCTz) + evA (T^g(w„ - LOq^UJu') - g{uJn,^n' - t^o)CT2 

z apii 



zc r 
'—A\dz^g{uJn - wo,a;„/) - V5(w„,w„/ -ujQ)dz 

-T^^\d{s, R; w„, ujk)g{s, R; ujk,0Jn') - g{s, R; ujn,ujk)a{s, R; Wfc, w„') 



(53) 



dA dg dg dA 



aR9p|| apii aR 

dd dq dc 



^E 



^ 2 



9A 9g dg dA 
9pil 9R 9R 9p|| 



dd 



do 



_9R9p|| 9p|| 9R 

where the quasiclassical matrix propagator is, as usualEalia 

dCp 



9p|| 9R 9R 9p|i 



g(s,R;w„,w„ 



G(p,R;w„,w„/) 



= 0, 



g -/ 

/^ 9 



and the order parameter is given by the self-consistency condition 



A(R) = giV(0)7r Y, J d'sJis, R; ^„ 



w„ 



(54) 



(55) 



Eqs. (|5^) and (|55| ) are the generalization of the standard quasiclassical theoryoEj to include terms giving rise to 
nonzero Hall conductivity. 

Before we linearize Eq. (^) and solve it to find the longitudinal and Hall conductivities, several comments should 
be made. First, the vector potential of the magnetic field enters the quasiclassical equation explicitly in contrast to 
the case of a normal metal (cf. Eq. (pl|)). This is readily understood if we notice that in the last term on the first 
line of Eq. (pSl) the matrix dzg — gdz has only off-diagonal elements, so that the term involving v^ only appears 
in equations for the anomalous propagator. It would seem that the second term involving the vector potential A 
(the last term in the second line of Eq. (53)) is present even in a normal metal since the matrix Uz'g + go'z has 



only diagonal elements, and, consequently, this term contributes only to the equations for the quasiparticle part of 
the matrix Green's function. However, in a normal metal in the presence of uniform electric and magnetic fields the 
response function is spatially uniform, and this term is irrelevant. In the superconducting state the spatial variation 
of the quasiparticle Green's function is due to the spatial dependence of the order parameter A in the vortex state, 
and this term describes the coupling of the current, induced by the spatial dependence of A(R), to the external field. 
Finally, the Lorentz force is accompanied by the matrix propagator in a combination azg + g^zi and therefore the 
Lorentz force does not act directly on the Cooper pairs. This result is perhaps not too suprising as in a reference 
frame associated with the center of mass the electrons have opposite momenta, and hence there is no net force acting 
on a pair. 



C. Linear Response 



We now use the approach given in Section [I C to linearize the basic equation in the external field. If we decompose 
the propagator 'g, self energy a and order parameter A into a leading order term go^ ag and Aq, and a part (denoted 
by index 1 ) linear in the applied electric field, the equation for the Green's function of a superconductor in a magnetic 
field reads 



"aAo 



go^z) + evA{azgQ - gao'z 



2(""^)a^(^^^° 



.90 0-2 



(o'ogo 

ie 
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9p|| 9R 9R 9p|| 
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9R apii 

go dao 



' Aq^o - ?oAo 



(56) 



= 0, 



while the equation for the response function g^^' is given by 



ivVg 
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9R 9p|| 
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da^ffg^ 
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9pil 9R ' 9R 9p|| 
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9R 9p|| 

9A(i) dgoi-) 

9pil dR ^ dK dp 



9pil 9R 

9q aA(i)' 



dgW aao(-) 



a^d) dgoi-) _^ 9go 5^(1) 
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9p|| 
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apii aR 



apii 9R aR ap| 



= 0, 



The rest of this work will be devoted to solving these two equations to determine the transverse electrical conductivity 
of a type-II superconductor in the vortex state. 

To calculate the response of a superconductor it will be convenient to modify the definition of current given in Eq. 
( p5| ) . Since the diagonal elements of the matrix propagator are related by Eq. (p^) , it is easy to check that the current 
can be written as 



j(R) = KeN{0) Y, I d'sv(s)(5i - 51), 



(58) 



where gi and gi are the diagonal elements of the response functixm-S^^^ In the standard quasiclassical approach the 
distribution function g also satisfies a "normalization condition"EJ't3 



i^k 



g{uJri,^k)gi^k,^n') = "(5, 



UJn.,UJ„> ' 



(59) 



In particular, using this condition for the leading order distribution function, which is diagonal in frequency (see Eq. 
(p7|)) we find 



dU^n) = -1- 



(60) 



However, it has to be emphasized that this normalization condition holds if and only if the gradient of the function 
g can be aajitten as a commutator of an operator with the distribution function, as is evident from the original 
derivationtj^EJ. It does not apply when terms responsible for the Hall effect are taken into account, since they have 
the form of an anticommutator of a matrix operator with the Green's function. Nevertheless, this normalization 
condition will prove useful in determining the quasiclassical Green's function of a superconductor in a high magnetic 
field at zeroth order. 
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IV. TYPE-II SUPERCONDUCTOR IN A HIGH MAGNETIC FIELD 

A. Model 

We consider a clean type-II superconductor in a magnetic field H close to the upper critical field iJc2- Again we 
consider a spherical Fermi surface, and impurity scattering is treated in the Born approximation. The condition for 
a superconductor to be in the clean regime is ^ ^ ^o- In fields not too far below the upper critical field the magnetic 
length A ~ ^Oj so that in the clean regime / 3> A. In type-II superconductors the spatial variations of the internal 
field become less pronounced as the superfluid density decreases with increased applied uniform magnetic field. As a 
result, near Hc2 internal fields can be assumed spatially uniform and equal to the applied figld and the vortex lattice 
can be modeled by an order parameter of the same form as the periodic Abrikosov solutior£3 

A(R) - ^ Cky'-y exp(-(x - A2fc^)2/2A2) = ^ Cky-y^„{x - K%), (61) 

h L- 

y y 

where $o(a^) is the lowest energy eigenfunction of the linearized Ginzburg-Landau Eq. (i.e. the eigenfunction of a 
harmonic oscillator with the Cooper pair mass M — 2m and frequency lOc). The vector potential of the magnetic field 
has been chosen in an asymmetric gauge ^(R) = {0,Hx,0). The periodicity of the coefficients Ck determines the 
type of the vortex lattice. Here, we do not consider a specific periodicity, the only assumption made is that there are 
flux lines in the system; this solution, therefore, can serve as a model for a rigid line liquid as well. 

B. Quasiclassical Equations in the Absence of an Electric Field and the BPT Approximation 

First we consider the leading order Eqs. (JS^) and neglect terms of order ujc- Then the elements are 

zvVg + A/t-AV=^A(/t)-^A*(/) (62) 

(2c^„ + v(V - 2ieA))f = 2iAg + -{f)g - -{g)f (63) 

T T 



(2c.„ - v(V + 2teA))f^ = 2iA*g + -{f)g - -{g)f\ (64) 



here angular brackets denote an average over the Fermi surface. The normalization condition, Eq. (p9|), can be used 
in this case, so that 

9^-.ff = -l (65) 

ff + ff = 0. (66) 

To solve these equations we employ the approach due to Bra.ndt, Pesch, and TewordtO, which was first used 
in the framework of the quasiclassical approximation by FeschE^O. In this method the diagonal elements g and 
g of the matrix propagator are approximated by their spatial averages, while the exact spatial form of A(R) is 
retained in determining the off diagonal functions / and / ' . The crucial observation is that the diagonal part of the 
Green's function is periodic in the center of mass coordinate R with the same periodicity as the order parameter. 
Performing apEourier decomposition of the full Green's function in the vectors K of the reciprocal fiux line lattice, 
these authorS showed that the Fourier components of the Green's function with K 7^ are exponentially small (by 
a factor exp(— A^iiT^)) compared to the component with K = 0. This component is, of course, the spatial average of 
the Green's function over a unit cell of the vortex lattice, which suggests the above approximation. 

The diagonal part of the distribution function depends on the amplitude of the order parameter, but not on its 
phase. The length scale for the suppression of the mean field order parameter amplitude by a single vortex is the 
coherence length ^Oi therefore near the upper critical field the order parameter is globally suppressed in the bulk 
of the superconductor. Consequently, spatial variations of the amplitude |Ap can be ignored for fields close to 
Hc2- On the other hand, as the phase of the order parameter changes by 27r around a single vortex, the rapid 
spatial variation of phase in the vortex state must be taken into account to determine the off diagonal elements 
of the quasiclassical propagator. After averaging over a single unit cell, the remaining spatial dependence of the 
amplitude |Ap is determined by the nonuniformity of the electromagnetic fields; the relevant length scale is the 
London penetration depth Al. Therefore, the BPT approximation works very well for superconductors in the London 
limit K = Xl/S,o 3> 1; even for materials with moderate values of k it remains valid over a wide field range below Hc2- 
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Numerical results obtained by BrandtES indicate that the BPT approximation works extremely well as long as the 
parameter (AA/w) < 0.3. Since the field dependence of the magnetic length is slow A « ^o(-f^c2/-ff)^ , this means 
that the approximation can be used over almost the entire region of linear magnetization, where the order parameter 
is suppressed. 

In all of the following g stands for the spatially averaged distribution function. To determine the functions 5, / and 
P we solve Eqs. (^3|) and ( |64| ) for the off diagonal elements of the matrix distribution function in terms of 5, and 
apply the spatially averaged normalization condition of Eq. ( pq ) to determine the diagonal part self-consistcntly. We 
introduce the impurity renormalized frequency 

U!ri=(^n + Tridi'^n)) (67) 

ZT 

and rewrite the equations for the off diagonal part of the distribution function as 

/ = (2(I3„ + v(V - 2teA)y\2tgA + ^(f)g) (68) 

/t = (2i5„ - v(V + 2ieA)y' {2igA* + -{P)g). (69) 



To proceed with this program we need to know the result of acting with the operator (2w„ ± v(V =F 2ieA)) on the 
order parameter. 

C. Operator Formalism 

Since the order parameter given in Eq. ( |6l[) is a superposition of the lowest energy eigenfuctions of a harmonic 
oscillator centered at different vortex cores, we introduce the raising and lowering operators 

a= ^[V^+ i{Vy - 2ieHx)] (70) 

v2 

a) = j=[^^-i{'sJy-2ieHx)]. (71) 

v2 

These operators obey the usual bosonic commutation relations [a, a^] = 1. We now interpret the Abrikosov solution 
as the ground state of this ensemble of oscillators A = |0). The higher eigenstates of the system are generated by the 
standard formula 



aV) = V^^H~T|n+l), (72) 

this operation excites oscillator states centered on each vortex line, so that 

I") = Y. Ck,e''''y'^n{x - A^ky) (73) 

ky 

Similarly we can introduce conjugate operators corresponding to A* — (0|, the raising and lowering operators for these 
states are now defined as 6 = (a)* and 6^ = (ct^)"^- Wide use of bosonic operators for the description of the vortex 
lattice has been hampered by the fact that, even though the wavefunctions corresponding to different oscillator states 
centered on the same vortex line are orthogonal, functions centered on different flux lines overlap, so thaL.different 
excited states as defined above are not orthogonal and the equations are non-local (see, for example, Ref.Ld). What 
makes this approach successful when combined with the BPT approximation is that this set of states is orthogonal 
in the sense of a spatial average 

' d^R(m\n)=A^6^^„, (74) 

where A is the spatial average of the order parameter. This condition is obeyed since the phase factor exp{ikyy) 
ensures that only functions centered on the same site contribute to the integral. Therefore if we are only concerned 
with spatial averages of physical quantities, the excited states of the order parameter can be treated as states of a 
harmonic oscillator. 
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To evaluate the result of acting with the gradient operator v(V — 2ieA) on the order parameter A we rewrite it in 
terms of the raising and lowering operators a and a^ 

v(V - 2teA) = ^^^ [ae-"^ - a^e"^] . (75) 

Then the result of the action of the operator (2a;„ + v(V — 2ieA}) on any mode \m) of the order parameter can be 
evaluated exactly. The technical details are given in Appendix |A|, here we give only the final result 

y — L OQ rrt 

(22„ + v(V - 2ie^))- Vi) = ^^^ y y 7?;^i™2e*('"^-'"i)'^|m + m2-mi), (76) 

?j Sin H ' ^ ^ ^ 



?7i2— mi— 



where 



(m — mi)!mi!TO2! v2 

2ia;„Asgn(a;„) 

Mn = ^—n > (78) 

wsmfc' 

W{u) = e""'erfc(-m), (79) 

and ly^™) is the rn-th derivative of the function W . Eq. (|77| ) is the main result of the operator formalism developed 
here; it allows further progress towards a solution of the quasiclassical equations to be made. 

D. Type-II Superconductor in High Magnetic Field 

Guided by the work of EilenbergerM and PeschES, we make an ansatz solving Eq. (Q) for an s-wave superconductor. 
This ansatz makes use of the fact that the term dependent on impurity scattering in the right hand side of the equation 
renormalizes the amplitude of the order parameter 

j = 2igD-^{Z)^){2Zjn + w{y -2ieA)\ A. (80) 

Since the order parameter A in this equation is the "ground sta te" o f the Abrikosov vortex lattice |0), the form of 
the function / can be obtained immediately from Eqs. (f6|) and (A1C|) 



/(s) = 2^g{s)D-\Z:r.)^^ £ -i=(-^)™e™^sgn(c^„))'"+'W^^'"n««)|m). (81) 

rn— 

Substituting this expression into Eq. ( pq ) , we find for the impurity renormalization of the order parameter 

I?(S„) = 1 - zV^-sgn(tj„) / d%(0; w„)VK(u„). (82) 

Using the corresponding Eq. (p3) we obtain /^(s) 



/t(.) = 2^g(^s)D-\Z^.^)^ Y. ^=(^)'"e-™*(sgn(^„))"+i W(™)(«„)(m|. (83) 

vsYnB ^ V"i \/2 

rn— ^ ^ 

Then we can use the normalization condition, Eq. (pq), to determine g (see Appendix H for details) 

5 = -isgn(Lj„)F(6',w„), (84) 

where 

r 2AA 9 n-i/2 

P(0,i„)- \-^^[——\'w'{u^) , (85) 
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and the sign has been chosen to give the correct expression in the normal state. Eqs. (^^ and (|8|) - (|8^) provide a 
complete self consistent solution of the quasiclassical equations for an s-wave superconductor in a nmgnetic field. A 
Green's function very similar to that given in Eqs. (r3) and (85) was obtained in the work of PeschEa by a different 



method. As in the microscopic theory, the order parameter is determined from the self-consistency condition given 
by Eq. (Es), which is in this case 



/^g7V(0)-V/ d0g(2„)i?-i(S3„)sgn 



K)W^K). (86) 

For the general case of finite mean free path and applied magnetic field a closed form solution of the self-consistent 
expressions cannot be easily found. However, with minor simplifications it is possible to obtain analytical results from 
this solution. Even though the dimensionless parameter (AA/w)^ in the Green's function given by Eq. (p3) is small 
in the region where the BPT approximation is valid, it appears with the weight (sin0)~^, so that a straightforward 
expansion is impossible. We will see, in fact, that the density of states is a non-analytic function of this parameter. 
However, while the full functional dependence of the Green's function on (AA/w)^ has to be retained, terms of 
higher order in this small quantity can be neglected in this functional form, provided that they do not result in more 
singular behaviour. Both the impurity renormalization of the frequency uJn and the renormalization of the order 
parameter depend on the weighted angular average of the Green's function g, which is non-singular as a function 
of the order parameter in the vortex state, this is related to the gapless character of the quasiparticle spectrum. 
Therefore, in determining the function P to leading order in (AA/u)^ the Green's function in the definition of 
impurity renormalization of the order parameter D (Eq. (p2)) can be replaced by its normal state value. Similarly, 
the renormalizcd frequency in the argument of the function W can be replaced by ujn +- sgn(a;„)/2r_-,The resulting 
expression for the renormalization function is identical to that obtained by Helfand and WerthamerE3. With these 
approximations, Eqs. (|8l|)-(|85|) describe a closed form solution. In the clean limit near the upper critical field of 
interest here expressions for the quasiclassical propagator can be simplified even further. Since in this regime I ^ A, 
and the renormalization of the order parameter is I? = 1 -I- 0(A//), to leading order Z? w 1. 

The anomalous Green's functions / and /^ are given as a Fourier series in the azimuthal angle </>, with the m-th 
component of the series coupling to the ?Ti-th excited state (or mode) of the order parameter A. Therefore in the 
presence of an external perturbation the mode with m = will couple to a scalar potential, the mode with to = 1 to 
a transverse potential etc. The function P given in Eq. (^) is related to the angular dependent density of states. If 
the Green's function is analytically continued into the upper half plane by letting iu;„ —> lu + i6 then the density of 
states 

N{lu, 9) = -N{Q)lmg{uj, 0) = iV(0)ReP(w, 6*) (87) 

is strongly angular dependent. For quasiparticles travelling parallel to the magnetic field, N{uj, 9) is gapped and 
BCS like, while in all other directions it is gapless. The total jdensity of states (DOS) Ns{uj), obtained by angular 
integration of the imaginary part of the function g, is gapless — |— while the residual density of states at the Fermi 
surface Ns{0) is a non-analytic function of the order parameterElIca 

/AA\2 /^/2v\ /AA 



, ,r /AA\-= /V^W\ /AAN^l , , 

iV..A.(0)[l-4(-) ln(^)+2(-)]. (88) 



The Green's function obtained here also reproduces the BCS Green's function if the limit iJ — > is taken, which 
suggests that it can be used to interpolate between the high-field and the low field regimes. We now have a closed 
form expression for the matrix propagator near the upper critical field up to the order (AA/v)'^, which we will use to 
determine the linear response of a superconductor to an electric field. 

V. LONGITUDINAL CONDUCTIVITY 

We begin by considering the longitudinal conductivity in the vortex state in the BPT approximation. We are 
concerned here with the transport coefhcicnts in the clean limit, and will neglect all contributions to conductivity 
of relative order (A//) compared to the most significant modifications upon entering the superconducting state. We 
again omit terms of order of cyclotron frequency. 
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A. The Response Function 

Since the electrical current given in Eq. ( |58| ) depends on ge — 9i—gi, we write the linearized quasiclassical Eq. (|57 
for the spatial average of this combination. Then the equations for the linear, in the applied electric field, averaged 
diagonal elements of the distribution function, and the equations for the anomalous functions are 

2evA(.g-.g(-)) 



ge = gi-gi = '^~ ^' " + i^^or\A*J - A*7H) + {iuo)-'{A^p - Ai/t(-)) (89) 

lUJQ 



+ (2ziior)-i i^iifhf (/!)/(-)) + ((/i)/t - (/i)/t(-)) 
~{2^^or)-' (^(7Rn - //(/(-)» + (Mn /i(/t(-))) 

2f2„ + v(V - 2ieA)\ h = ievA(/ + /(-)) + iA(gi - 51) + iAi(.g + g{-)) (90) 

+*(2r)-i((A>(.g + .g(-)) - {f)-g, + (/(-))gi) 

2r2„ - v(V + 2ieA)\ fl = ievA(/^ + /^(-)) + zA*(gi - gi) + iAl{g + g{-)) (91) 

+^{2r)-'[{fl){g + g{-)) (/t(-)).gi + (p)g,). 

The notation used here is identical to that of the previous Section, and the frequency il„ is defined as r2„ = u;„ +<!;—. It 
is possible to identify the different contributions to the right hand side of Eq. (p9h . The first term is the quasiparticle 
contribution to the current, this term determines the response function in the normal state, and, with the modified 
Green's function, describes the contribution of quasiparticles to the current in superconductors. The other terms on 
the right hand side exist only in the superconducting state. The first two of these involve the modification of the 
order parameter A, and can be associated with the motion of the vortex lattice under the influence of the applied 
electric field. The remaining terms mix the contributions of the quasiparticles and the Cooper pairs. It will be shown 
below that the most relevant contribution from these terms is due to the additional scattering of the quasiparticles by 
dynamicaLfluctuations of the order parameter, similar to the processes described in the dirty limit by the Thompson 
diagramsO 

The quasiparticle contribution to the response function gi—gi can be determined immediately since the unperturbed 
functions g and g{—) arc known from Eq. (|84|). To evaluate the other contributions to the response function, Eqs. 
( P0| ) and (|9l| ) have to be solved for Ai and (/i), as well as for the conjugate quantities A* and (/|). As before, here 
we determine the complete functional dependence of the response function on the order parameter to order A^ , and 
neglect corrections that vanish faster than this as A decreases. Since both Ai and /^ can be expanded in a complete 
set of functions \m) and (to|, which are normalized by A^, see Eq. ([7^), it is sufficient to determine the expansion 
coefficients to zeroth order in A. Therefore in Eqs. ( po|) and (^ij) we can replace the functions g, gi and gi by their 



normal state values. With these simplifications Eqs. (|89D-(pi[) can be solved exphcitly for Ai, /i and the "daggered" 
functions. 

B. Quasiparticle Contribution 

Two different effects modify the quasiparticle contribution to the current relative to the current in a normal 
metal. First, the difference g — g{—) is modified relative to its normal state form, and, second, as the impurity 
renormalization of the frequency ujq depends on the unperturbed Green's function it is also affected by the opening 
of the superconducting gap below the upper critical field. In the normal state, the difference g — g(— ) vanishes in 
the outside frequency region, for a type-II superconductor this difference is of order A^. Therefore in a calculation to 
lowest order in A^ the renormalized frequency can be replaced by the bare frequency in this frequency range. On the 
other hand, in the intermediate frequency range, where the difference of the unperturbed Green's functions 5 — .g(— ) is 
of order 1, it is important to keep the full dependence of the renormalized frequency on the order parameter. Further, 
as the contribution from the outside region is proportional to to, but not r, it is of order (AA/u)^(A//) and negligible 
compared to the contribution from the intermediate frequency range. This situation is not unusual when comparing 
different contributions to the conductivity. Two dimensionless quantities involving the frequency of the external 
electric field appear in our analysis. The first, wr, usually comes from renormalization of the bosonic frequency ujq 
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in the intermediate frequency range. The second, (A(D/u), appears when the response functions are expanded in the 
external frequency since the argument m„ of these functions involves the frequency in the combination (Auj/v), see 
Eq. d?^). In the dc response only terms linear in uj contribute to the absorptive part of the conductivity. Therefore, 
as the ratio of the two dimensionless parameters is of order (A/Z), we keep terms of order lot while neglecting those 
of order (Aoj/v). In the quasiparticle contribution then the only relevant terms arise from the intermediate frequency 
range. 

Since the quasiparticle spectrum is gapless in the high field regime, the response function varies slowly over the 
scale Lu ^ T, and the frequency sums can be evaluated easily, see Appendix ^ We find the quasiparticle contribution 
to the current 



1 



Jgp 



iv(o)e^w^A / shrede 



(P - 1) + iuiT F + ((1 - P)) - Vt^W" 



A 



2AA 



I sin 6 J \v sin 



A 



pj 



(92) 



here all the functions are evaluated at w = 0. For w = the argument of the function W' in Eq. ( p5|) is purely 
imaginary, and the function P is purely real. It follows that the first term in Eq. (92) contributes to the non- 
absorptive part of the conductivity; it is the remnant of the Meissner effect in a type-II superconductor in a magnetic 
field. The remaining terms contribute to the absorptive part, and the transport current can be written as 



hp = ]N{0)eWi: f sin3 
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The first term in Eq. ( P3| ) is the direct modification of the quasiparticle current on entering the superconducting state 



iqpi = ]N{0)e^v^Tl] I sin 
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where the normal state conductivity (t„ was defined in Section [ID . The contribution of small angles sm6 < (A/l) to 
the angular integrals is of higher order in (A/Z) and can be neglected. For larger angles the argument of the function 
W can be set to zero since A/ 1 <C 1. Then the integration is easily carried out, expanding the resulting elliptic 
integrals for small values of the parameter (AA/ii), we find that the correction to the conductivity from this term 
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is negative. In the superconducting state in addition to the scattering of quasiparticles by impurities, quasiparticles 
are scattered by the vortex lattice. At a vortex core a quasiparticle can undergo Andreev scattering into a hole and 
a Cooper pair with no energy cost. This additional scattering process reduces the quasiparticle contribution to the 
current. 

The second term in Eq. (p3) arises from renormalization of the scattering time r in the vortex state. It can be 
written as 



iqp2 
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E = iiV(0)e2«Ve//E, 
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where the scattering rate 
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The quantity N{0){P) evaluated at w = is the residual density of states in a superconductor N^, see Eq. (pq). Hence 
this term describes the effect of the change in the density of states on the scattering rate of the quasiparticles. Below 
the transition, as the superconducting gap opens, the residual density of states at the Fermi surface is suppressed 
compared to the density of states in the normal state; consequently, the effective scattering rate is smaller and the 
effective mean free path is larger. The angular integral of the function P can be evaluated to leading order in {A/l) 
and expanded in (AA/w) in similar fashion to the integral analyzed above, we obtain, in agreement with the result of 
Eq. (pi), the effective scattering time 
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and the contribution to the longitudinal conductivity 
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Since (AA/f) <IC 1, the logarithmic term dominates near the transition and this contribution is enhanced relative to 
the normal state value. The last term in Eq. (BS) contributes at order (AA/u)^(A/Z). 



C. Dynamical Fluctuations of the Order Parameter 

To compute the contribution of all the other terms in Eq. ( p9| ) to the current, we have to solve Eqs. ( pO| ) and ( ^l| ) for 
the linear, in the electric field, correction to the order parameter and determine the functions /i and /|. As discussed 
above, the functions gi and gi can be replaced by their normal state values to the order to which we work. To use the 

operator formalism we need to evaluate the effect of acting with the differential operator, 251„ + v(V — 2ieA) , 
on the unperturbed function /. In the clean limit 



/ = 2z5 



2w„ + v(V-2ieyl) 



Then the two differential operators can be separated 
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and Eq. (90) becomes 
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Since the self consistency condition requires that 

Ai = ^giv(o)r^(/i) 



(103) 



and Eq. (102) is a standard Fredholm type integral equation, it is clear that, since the function / given in Eq. (^ 



is a Fourier series in 4>, the angular average of the product / cos 4> projects out only the component proportional to 
exp(«(^), the first excited mode of the order parameter \m = 1). This implies that the linear, in the electric field, 
change in the order p aram eter involves only the first excited state, and has the form Ai ~ C |l), w here C is to be 
determined from Eq. (102), sipularly, A^ — C{\\. This result, which was anticipated in Section IVD , is in agreement 
with that of Caroli and MakiE3. The contribution to the response due to the dynamical fluctuations of the order 
parameter is given by the second term in Eq. (|89|). Using the functions / and /^ from Eqs. ( po|) and ( |9l| ) and the 
orthogonality condition given in Eq. ([74|) , we find that 
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and the contribution to the longitudinal and transverse electrical current is 
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We see that the "odd" part of the dynamical fluctuations of the order parameter gives rise to a contribution to 
the longitudinal resistivity, while the "even" part contributes to the Hall current. Both of these contributions are 
proportional to the "vertex" function 



y(c^o) = T 



gW'-g{-)W'{-) 



IWo 
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describ ing t he coupling between the electric field and the excited mode of the order parameter. 

Eq. (102) is solved in Appendix ^, we find that the last two terms result in small, in (A/Z), contributions, and the 
amplitude of the fluctuations C is given by 



C = ^evA 
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The denominator on the right hand side of Eq. (109) is the propagator of the first excited mode of the order parameter. 
In general, the zeroes of this propagator correspond to the spectrum of propagating modes of the order parameter. 
In our case the transverse perturbation due to the vector potential of the electric field couples to the first excited 
mode of the order parameter, which is damped, i.e. there is a finite energy gap in the spectrum of these excitationSpajt 
zero frequency. The response to a scalar potential is quite different, there is a propagating mode at zero frequencju. 
Since the dynamical fluctuations of the order parameter are driven by the electric field, the coupling to the excited 
mode in the numerator of Eq. (|109D is also proportional to the vertex function defined in Eq. (108). 
Evaluating the sums in Eq. (|109|) we find (Eq. (|C6|)) 



C = -C = 
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and therefore there is no contribution to the transverse current due to the fiuctuation term, as expected, 
contribution to the longitudinal current can be evaluated from Eq. (106), it is 
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and the contribution of the dynamical fluctuations of the order parameter to the longitudinal conductivity is 
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D. Thompson Contribution 

We now consider the remaining terms in Eq. (p9). To evaluate their contribution to the longitudinal current we 
need the explicit expressions for the angular averages of the unperturbed anomalous Green's function / and /^, and 
the linear, in the electric field, corrections /i and // to the distribution function, these are obtained from Eq. (102) 
and its daggered counterpart. 

Only one of these terms, the term involving the angular average of the functions /i and f\^ gives a contribution to 
the conductivity at the order considered here. A typical term is given by 
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where the functions under the integral depend on the angle 6' , and 
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2zM7„sgn(17„) 
V sinfl 



(114) 



in analogy to Eq. ([78|). As there is an additional factor of the scattering time in the denominator, it might be 
expected that this contribution is small. However, in the intermediate region the renormalized frequency wq ~ 1/''" 
and, since g and g{—) have opposite signs, the contribution of gW' — 9{—)W{—) is of order one. Therefore the first 
term contributes to the conductivity at the same order as the corrections found previously. On the other hand, as 
g + g{—) = in the intermediate region, the second term does not contribute to the current. In the outside region 
both terms give contributions to order {K/l) which can be neglected. The contribution to the current from Eq. (|113|) 



and the corresponding term involving (//)/ (which is obtained from equation (113) by replacing e '"^ with e*"^) is 



j^,i = !l!7V(0)eVA(— )'-^i V (\inede {w' + w'{~)] 



The Thompson-like contribution to the conductivity is given by 
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In his original work ThompsorO found that there is a contribution to the conductivity in the dirty (/ <C ^o) limit 
due to scattering of quasiparticles by the dynamical fluctuations of the order parameter. The main contribution in 
the dirty limit arose from the outside region; the contribution of the intermediate region was smaller by a factor 
{1/^q). The result obtained here for the clean limit is consistent with this picture. The term contributing to leading 
order is proportional to the angular average of /i and exists only in the presence of the excited mode of the order 
parameter as it depends on the angular average of /i. As expected when (Z/^o) ^ Ij the relevant contribution comes 
from the intermediate region. In the presence of a transport current the vortex lattice moves, and individual vortices 
are deformed. As a result, additional scattering of quasiparticles by the vortices gives rise to a negative contribution 
to the conductivity given in Eq. (116). 



E. Longitudinal Conductivity 

The longitudinal conductivity of a clean type-II superconductor in the mixed state is obtained by combining the 
results for the quasiparticle current from Eqs. (|95| ) and (|9^), the current due to the dy nam ical fluctuations of the order 
parameter from Eq. ( |112| ) and the current due to the Thompson terms from Eq. (116). We notice that reduction 
in the quasiparticle contribution to the conductivity due to additional scattering off the ground state of the vortex 
lattice (Eq. (pa)) and the excited modes of the order parameter (Thompson terms) is compensated to order (AA/w)^ 
by the increase in the current due to dynamical fluctuations of the order parameter The conductivity then is given 
by Eq. (M) 
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that is the modification of the longitudinal conductivity upon entering the superconducting state is determined solely 
by the increase in the effective mean free path due to the suppression of the density of states at the Fermi level as the 
superconducting gap opens. The increase in the mean free path is a non-analytic function of the order parameter. 



VI. HALL EFFECT 



A. Stability of the Leading Order Solution 



In determining the density of states and the longitudinal conductivity we have neglected terms of the order of 
cyclotron frequency not only in the linearized quasiclassical equations, but also in the leading order equations (63)- 
(p4) . To investigate the behavior of the transverse conductivity the gradient terms in Eq. (pq) have to be taken into 
account, and the solution for the propagator g at zeroth order in the electric field has to be obtained to order lUc- 
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Instead of attempting to solve equation ( pq ) in full, we will show here that the solution obtained in Section IV is still 
valid when terms of order of the cyclotron frequency are included in the equations. 



We saw in Section III that, as the matrix combination azg + gdz has only diagonal elements, the Lorentz force 
acts only on the quasiparticle (diagonal) part of the propagator. The function g given in Eq. ( p4| ) does not depend 
on the azimuthal angle (p, and, therefore, there is no correction to this function from the Lorentz force term. Next 
we observe that the term involving the gradient of the propagator in the third line of Eq. (|56| ) is proportional to 
the same combination of matrices as the Lorentz term. Since in the BPT approximation the function g is replaced 
by its spatial average, this term vanishes. For an s-wave superconductor the order parameter A is constant at any 
point at the Fermi surface, and its derivative with respect to the components of momentum parallel to the Fermi 
surface vanishes, which means that the last term in the fourth line of equation (|56| ) can be ignored. The momentum 
derivative of the self energy due to impurity scattering vanishes for the same reason. 

We now consider the remaining terms in Eq. (pq). Omitting the subscript since in this Section we only consider 
functions at leading order, we write the first of these terms in the matrix form 



M= - 
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The off diagonal elements of this matrix are proportionl to the trace of the quasiclassical propagator and vanish in 
accordance with the normalization condition. The contribution from the term M to the equation for the quasiparticle 
part of the distribution function is 



Mil = - 
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dA dp dA* df 



dm 9p|| dR 9p|| 



(119) 



To spatially average this term and determine its contribution to the diagonal part of the propagator we need to recast 
the gradient operators in terms of raising and lowering operators a and a^ and the azimuthal and polar angles (j) and 
9. We find 
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Here the hat denotes the gauge invariant gradient 
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and the operator with the plus si gn a cts on A* while the operator with the minus sign acts on A. A direct check using 
the solution obtained in Section IV shows that the terms breaking gauge invariance vanish after spatial averaging, 
consequently, the gradient can be replaced by its gauge invariant counterpart, as expected for an operator acting on 
the order parameter. For the ground state of the vortex lattice a direct check shows that this term does not result in 
any correction to the unperturbed propagator. The contribution from the term involving the spatial derivative of the 
self energy vanishes in complete analogy to the term just discussed as their structure is identical. 

Therefore the solution of the quasiclassical equations obtained in Section IV also satisfies the quasiclassical equations 
when terms of the order of cyclotron frequency are taken into account. 



B. Linearized Equations for the Transverse Response 

We now consider the linearized Eq. (p7[). In the regime when lOcT ^ 1 terms of the order of cyclotron frequency 
can be included in the calculation of the response function perturbatively. We therefore solve for the linear, in the 
cyclotron frequency, corrections to the averaged response function ge = <7i — .91 obtained in the calculation of the 
longitudinal conductivity in the preceding section. 

Since the Hall conductivity in the normal state is proportional to the square of the scattering time r, we can expect 
that the most relevant contributions to the transverse conductivity in the vortex state are also proportional to t^, 
other contributions to the Hall effect are smaller by a factor (A/^). Therefore we keep in the equations only terms 
that contribute to this order. If now 5g is the part of the propagator linear in the cyclotron frequency, we arrive at 
the following equation for the function 5ge = Sgi — 6gi 
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Here (71 — gi is given by Eq. (B9[) . There are now two distinct contributions both to the term involving the fluctuation 
of the order parameter and to the Thompson term. One reason these terms contribute to the transverse response 
is that they give rise to additional scattering due to dynamical fluctuations of the order parameter induced by the 
electric field, as we saw in the previous section. When the quasiparticle trajectories are bent by the magneti c fie ld, this 
additional scattering renormalizes the Hall conductivity. This effect is contained in the first term in Eq. ( |l22| ), since 
the function gi — gi contains the contributions of the Thompson terms and the fluctuations of the order parameter 
induced by the electric field. The other contribution to these terms is due to fluctuations of the order parameter 
induced by the Lorentz force, these fluctuatio ns r esult in corrections to the transverse conductivity and are contained 



in the terms involving SAi and Sfi in Eq. (122). Finally, the terms involving the gradient of the order parameter 
contribute to order r^ for the same reason that the Thompson term contributes to the longitudinal conductivity, 
namely, that there is an additional factor of the scattering time r in the amplitude of the order parameter fluctuations 
C and in the functions /i and fl, so that the overall contribution is of order t^. This anomalous contribution to 
the transverse conductivity arises because the gradients of the order parameter created by the moving and deformed 
vortex lattice act as driving forces (analogous to the Magnus force) in the transport-like equations. The remaining 

terms in Eq. (M) contribute at higher order in I A/ 1 j . 

The equations for the corrections to the off diagonal elements of the matrix distribution function are 

[2a. + v/(V - 2ieA)]df = -^a(v - 2ieA) (/ - /(-)) + lASg, + iSAi(^g + g(-)) (123) 

i dA d f _ \ i 9Ai 9 / \ 

+ 2aRa^l^^+-^V + 2 9R9^l^ + ^(^J 

and 

[2(7„ - vj(V + 2ieA)]Sp = -^^(v + 2ieA) (/^ - /t(-)) + iA*dg, + iSA^ (.g + .g(-)) (124) 
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In the normal state there is no angular dependence to the unperturbed funct ion g , and also 51+51 = 0, so that the terms 



in the last line of each equation vanish. We can now solve Eqs. ( |l23| ) and (124) to determine the contributions to the 
dynamical fluctuations of the order parameter induced by the Lorentz force. We can then evaluate the contributions 
to the transverse conductivity term by term. 

C. Hall Conductivity 

The quasiparticle part of the response function is 

Sg^P = -2evAsmesin^'-^ ^^ ~ j^~^^ (125) 

iujQ iujQ 

The contribution to the conductivity from the intermediate frequency range is readily evaluated; the correction to the 
transverse conductivity due to additional scattering off the vortex lattice 

Aalf = -6a„c^,T(AA/«)2 (126) 

and the contribution to the Hall conductivity due to the modification of the scattering time 
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(127) 
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In addition there is a quasiparticle contribution to the Hall conductivity from the outer frequency range which is 
formally divergent 



jr' = -^7V(0)e^2A^ / sin^' 
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Since P' evaluated at zero frequency is purely imaginary, the second term describes a small correction to the Meissner- 
like term. The first term in this equation, on the other hand, has no physical meaning and must disappear from the 
final expression for the current. 

It is in fact cancelled by the contribution of the fluctuations of the order parameter 

4 ^2^.eu.CN{Q)^-— d s.—^T}^ ^r^, (129) 

from the outer range, where the first term in the expansion of the vertex function 

out ^ "' ui„>Q 

is formally divergent. The remaining contribution from the outer range is obtained by expanding the coefficients C 
and C in small quantity ujt, it is 

jf = 6a„(c.,T)(^)'£;. (131) 



This contribution is cancelled by that of the intermediate frequency range in Eq. (129), so that there is no net 
contribution to the transverse conductivity due to the dynamical fluctuations of the order parameter driven by the 
electric field. This result is consistent with the predictions of time dependent Ginzburg-Landau theorylil. The terms 
considered so far in this Section correspond directly to those contained in TDGL, which is an effective theory treating 
only the fluctuations of the order parameter, while the quasiparticle contribution is taken to be at the normal state 
value. In the TDGL approach the Lorentz force has no effect on the dynamics of the order parameter, and there is 
no correction 6Ai due to this force. 

In the present analysis, however, the equations for the quasiparticle propagators and the amplitude of the order 
parameter fluctuations are coupled, so that even though the Lorentz force does not appear explicitly in the equation 
for Sf, it introduces changes in the diagonal part of the distribution function gi — gi and therefore brin gs ab out 



further modification 6Ai of the order parameter. To find this contribution we have to solve Eqs. ( |123| ) and (124) for 
the changes in the order parameter SAi — SC\1) and SA* — SC{1\. The solution, given in Appendix O, follows the 
same steps as in the calculation of the longitudinal conductivity. We find that only the term involving 6ge contributes 
at the order to which we work, and 

5C = SC = ieAAV2{u!T){ujcT). (132) 



and, since we saw in Eq. (107) that the "even" part of dynamical fiuctuations of the order parameter contributes to 



the transverse part of the conductivity 



jfl = Qan{uJcT)(^y. (133) 



The dynamical fluctuations of the order parameter driven by the Lorentz force tend to increase the transverse con- 
ductivity. 

Similarly, there are two parts to the Thompson terms: one is due to the longitudinal response, while the other is 
due to the linear in the Lorentz force corrections to the off diagonal distribution functions. The contribution to the 
transverse conductivity due to the first of the Thompson terms is the longitudinal contribution multiplied by lucT, 

a™ = -6a„(c.,r)(^)'. (134) 

To determine the contribution to the current from the second Thompson term, we use the angular averages (df) and 
((5/^) given in Appendix O, to find that its contribution doubles that given in Eq. (134), so that the total contribution 
of the Thompson terms to the transverse conductivity is 
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In addition to the scattering of quasiparticles by the fluctuations of the order parameter induced by the appUed 
electric field, which tends to reduce the current, quasiparticles also undergo additional scattering off the dynamical 
fluctuations driven by the Lorentz force, which again tends to reduce the transverse response. In this sense, the 
Lorentz force results in anisotropic scattering of the quasiparticles by the fluctuations of the order parameter. 

The l ast tw o terms, the gradient terms in Eq. ( |122| ) have a structure identical to that of the terms considred in 
Section VIA. Their contribution to the current can be evaluated using the operator approach as shown in Appendix 
W. First, consider the term involving the gradient of the dynamical fluctuations of the order parameter Ai and A*. 
These fluctuations involve the first excited mod e of the order parameter, which, when they are acted upon by the 
annihilation and creation operators in Eq. ( A17), gives terms proportional to the ground state and the second excited 
state of the order parameter, respectively. Then spatial averaging projects out the same modes from the functions / 
and /^ given in Eqs. ( pi|) and (p3|). In the second of these terms the gradient operator acts on the ground state of 
the vortex lattice. A, which the creation operator promotes to the first excited state. Spatial averaging now projects 
out the first excited component from the functions /i and fl . The contribution to the transverse current due to these 
terms, is found to be 
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the contribution to the transverse conductivity due to these terms is given by 



< = 12a„o;.r(^)^ 
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The induced gradients of the order parameter enhance the transverse conductivity. 

The total transverse conductivity is the sum of all the contributions considered here. Wc find that t he m odification 
of the quasiparticle Hall current due to additional scattering off the vortex lattice given in Eq. (126) is exactly 
compensated by the enhancement of the transverse current due to Lorentz force driven fluctuations of the order 
parameter obtained in Eq. (133) Th e Th ompson contribution due to additional scattering by the deformed and 
moving vortex lattice is given in Eq. (135) and is cancelled by the enhancement of the transverse conductivity due 
to the forces generated by gradient of the excited mode of the order parameter found in Eq. (138). As a result, the 
behavior of the transverse conductivity axy is determined solely by the modification of the effective elastic scattering 
time Teff and is given by equation (127) 
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For the dc conductivity this change is due to the decrease in the number of states at the Fermi surface available for 
scattering as the superconducting gap opens. 

VII. CONCLUSIONS AND DISCUSSION 

We now plot qualitatively the longitudinal resistivity (Figure IT]), the transverse conductivity (Figure 0) and the 
Hall angle (Figure ra) as functions of the applied magnetic field for Niobium. The order parameter, whijch is linear in 
the applied magnetic field in the high field regime, is given by the expression due to Maki and TsuzukiE3 
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and the values of the superconducting material parameters were taken from Refs.EErEfl. The longitudinal resistivity in 
Figure ^ has a pronounced increase in slope as a function of the magnetic field below the superconducting transition 
due to the logarithmic dependence in equation (117). The transverse conductivity shown in Figure || is enhanced 
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below the upper critical field and has negative curvature in the high field region. The negative curvature arises from 
the competition between the enhancement due to the increase in the effective mean free path and the linear decrease 
of the cyclotron frequency with the field; the Hall conductivity is substantially enhanced when compared to the linear 
decrease expected from the normal state behavior. While the transverse conductivity is proportional to the square of 
the scattering time, the Hall angle 

tandn = a-xy/(^xx = ^cTe// (141) 

is only linearly dependent on the scattering time and the corresponding nonlinear dependence on magnetic field is 
weaker, as can be seen in Figure |3[ Finally, as the transverse resistivity 



P^y ^ ~ — 7^"^~ ~ ^xy/olx (142) 
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is independent of the effective scattering time, it remains linear in magnetic field upon entering the superconducting 
state jacith the same slope as in the normal metal. This behavior is to be contrasted with that of Bardeen-Stephen 
modeH, where the resistivity is modified and is linear in the magnetic field, but the Hall angle obeys the same linear 
law as in the normal state. The Nozieres-Vinen theor};^, on the other hand, which predicts that the Hall angle should 
be constant in the flux-flow regime below Hc2 at variance with the result of this work, also finds that the transverse 
resistivity is identical to that of the normal state, although the individual components of the conductivity tensor are 
quite different from those found here. |-| 

A comparison can be made with the experimental data of Fiory and SerinQ on high purity Nb. These experiments 
find a transverse resistivity in the flux-flow regime which is linear in the applied magnetic field over a wide range of 
fields below Hc2- The Hall angle, however, flattens or even increases above its value at Hc2 before decreasing at lower 
fields. These results are more suggestive of the behavior given here than the original interpretation given in terms of 
the Nozieres-Vinen theory. Also, the longitudinal resistivity found in Rcf.Q has a distinct increase in slope just below 
the upper critical field, which is consistent with the behavior discussed above. Detailed comparisons with the results 
of this work are difhcult to make, since the authors of Ref.B used a high current density to reduce the pinning effects 
and achieve the flux flow regime; as a result, the magnetoresistance is significant and the longitudinal resistivity in 
the normal state varies with magnetic field. We find the qualitative agreement with the experiment encouraging and 
suggest that more experimental work is needed to make a more detailed comparison with the theory. To conclude, we 
have presented here a new approach to the calculation of the transport coefficients of a clean type-H superconductor 
in the vortex state in the high-field regime and used it to determine the Hall conductivity and the Hall angle of an 
s-wave superconductor in this regime. We find that the field dependence of the Hall conductivity in the high field 
regime, which is non-analytic, is entirely due to the change in the density of quasiparticle states at the Fermi level in 
the superconducting state. At the same time we find that the field dependence of the transverse resistivity below the 
upper critical field remains unchanged. 
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APPENDIX A: THE OPERATOR FORMALISM 

If \m) is the 771-th excited mode of the order parameter we have 

[2a;„ + vj;^(V-27e^)]"^|?77) =sgn(w„) / expf -(2(I;„ + vf(V - 27e^))sgn(w„)f Jdt |?77) (Al) 

= sgn(a;„) / exp(-2a;„sgn(a;„)i) expf -=— [ae"*'^ - a'fe"^]sgn(u;„)i jdi |?77). 



We use the operator identity 



e^+^=e-4e^e-^[^'^l, (A2) 
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where [A, B] = AB — BA denotes a commutator, to separate the creation and annihilation operators and rewrite Eq. 
&) as 



[2w„ + vf(V - 2ieA)]~'^\m) = sgn(, 



uin) I diexp(-2w„sgn(w„)i -— 5— 

4A^ 



X exp 



vt sin 



v^A 



sgn(a;„)e"^a^' 



exp 



wising 

V2A 



sgn(a;„)e *'^a |to 



We now write the exponentials as infinite series in powers of the arguments to find 

°° ™ foo ,,2o;,i2, 



Z""" u^sin (-1)™! 

[2ii„+VF(V-2ze^)]-V> = = I" I" / dtexp(-2<i„sgnK)* 77^—^') , , e'^" 

•■^ -^^ Jo 4A"^ mi\m2'. 



rn2— rni— 



iitsin^ 



V2A 



?-nl+?7T.2 



x(sgn(tj„)) 
In the integral, the parameter t can be replaced with a differential operator 

t 
and the integral can be evaluated 



t\"^2 / \m 



(aT) ^(ani"^)- 
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.,2 „:„2 , 



/ diexp — 2ijj„sgn((jj„)i r-r^ — t 

Jo L 



77il+m2 



1 a 



wsin0 

V2A 



4A2 

mi+m2 /.oo 



ml-\-ni2 



vtsinO 

V2A 

(iiexp(— 2a;„sgn(a;„)i 



w^ sin^ I 
4A2 



■t') 



■ysin^^ y/2 
where W{u) = e~" erfc(— iu), H/(™) denotes the m-th derivative and 

_ 2ia;„Asgn(w„) 



The main result is 



V^A 



V sin 6 



CX3 m 



[2W„ +Vf{\/ - 2ieA)]-^\m) = ^^ ^ ^ £,mim2g^(m2-r„l)0|^ + ,„2 - mi). 



V sm 



7712—0 77^1—0 



where 



^7771,^2 ^ ^iz!\/i^l™l^^(„l)™l(_ * )™l+™2(sgn(^,^))"l+™^ + V(™l+™^)K). 



(m — miy.milm2l 



V2' 



We make extensive use of two special cases of Eqn.( A8): 



(215,7 + v(V - 2*e^))-i|0) = ^^ V (-^)'»e-'A(sgnK))"+iW^(")K)k 

wsm6'f-' V^ V2 



777 = 



and 



(22}77 + v(V ~ 2^eA))-^\l) = ^ ^ (_^)™(sgnK))'"+' 

\/^^me""'^H/('")(M,7)|TO + 1> + (^)(sgnK))e*("-i)^M^("+i)K)|m) 

v2 



(A3) 



(A4) 



(A5) 



(A6) 



(A7) 



(A8) 



(A9) 



(AlO) 



(All) 
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Equations for the daggered quantities are obtained by replacing the phase imcj) by its conjugate —imcl), changing the 
sign of {i/^/2), and using a bra vector instead of the ket vector. 



To determine the quasiclassical Green's function g we need the spatial average //t. Using Eq. (A4) we have 



//t = / d^Rfp = -Va2 / dhdt2 exp -2w„sgnK)(ti + ia) - 



v"^ sin^ ( 



4A2 



■(^1+^2)' 



«v^5^ 



2AA 



W 



,/2iLj„Asgn(cj„) 



(A12) 
(A13) 



Eq. (|84| ) obviously follows from the last line. 

In the calculation of the transverse conductivity we will need to rewrite the gradient operators in terms of creation 
and annihilation operators. Since a gauge invariant gradient can be written as 



d__ J_ 

dx ~ V2A 
d I 



1 



By V2A 



a + a^ 



V2A 
i 



V2A 



6 + &^ 



and the the momentum gradient in the direction parallel to the Fermi surface is 



-ee- 



1 ^ d 



9p|| p dB psin 
where eg and e<j, are the unit vectors in d and </> direction in the spherical coordinates, we obtain 



(A14) 
(A15) 

(A16) 



d_ d _ 1 

5R9p7 ~ V2pA 
1 



('ae-*"^-a1'e*'^]cos( 



V2pA 



he"^ - h'^e-"' 



d i 

I 

89 sin 
.A 

ae 



' V /I 



sm 0\ / ( 



(A17) 



APPENDIX B: FREQUENCY SUMS 



The sum of the values of a response function at Matsubara frequencies iujn = (2»t. + l)TriT in the upper half plane 
can be written as an integral 



^E^(-") = i-y^tanh( )kH. 



(Bl) 



If the response function varies slowly over the scale to ^ T, and the tangent can be replaced with a step function, so 
that 



Ty^K{iuj„)K — :flimF(w)+ lim F{uj) - 2F{0)] , 



(B2) 



n=0 



where K^ui) — (dF (uj) / duj) . 

First consider the sum that appears in the quasiclassical contribution to the longitudinal current 



^-E 



{g-gi-)) 



IWO 



(B3) 



Since the frequency Wq can be replaced with bare frequency in the outer frequency range, but is renormalized in the 
intermediate range, we consider the sum separately in the two regions. In the outside region, transforming the sum 
in the lower half plane into a sum over the frequencies in the upper half plane. 
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s, 






2 

UJQ 



w„>0 w„>0 

after analytic continuation and expansion in w. Using Eq. (p32|) we obtain 



9 -^ 

TT 



(P-1) 



Ad- \dP 



V sin / did 



(B4) 



(B5) 



where the values of the functions are computed at w = 0. To evaluate the sum in the intermediate frequency range 
to leading functional order in A^, we write 



wo P + P( — ) '^° 1 

^"* ^ ~'^,^^o + mT)i{p) + {p{-))) = '^'^/ij;^TW^^wi+lPH))' 



(B6) 



Adding and subtracting the contribution of a normal metal, so that the remaining sums are convergent at high 
frequency, we obtain after analytic continuation 



^it 



(--') 



l + ( 1-P 



(B7) 



and 



5 = 



iA \ / 2AA \ 2 / A 



(p-i) + -r[p + ((i-P))-v^^"(-^)(^)7-^) 

L V i sm 6' / V i; sm t' / V / sm 6^ / 



P' 



(B8) 



The vertex appearing in the calculation of the dynamical fluctuations of the order parameter is proportional to 

gW'-g{-)W'{-) 



E 



lUjQ 



(B9) 



Since the amplitude of the dynamical fluctuations only has to be evaluated to zeroth order in the superconducting 
order parameter, it is sufficient here to replace the Green's function in the renormalized frequency by its normal state 
value. Since this sum is well-behaved at high frequency, we easily obtain 



^, _ y^ gW - 9{-)W'{~) _ ^. 



\lUJn + I T ILOnJ ^ — ' V 



t^„>0 



After analytic continuation we find 

_ -2i 2A 



Aw 



1 — iujT V sin 
1 1 



S *•"'»•■' + (^)'^'"<»") 



TT 1 — iojT 



W 



w„>0 



I sin 9 / V sin 



Au; ^„/ zA 



I sin 6 



(BIO) 



(Bll) 



The sum in the fluctuation propagator in Eq. ( |10S| ) is easily evaluated in a similar fashion 

1, 



J2 [*5SgnK)W^(u„) - ^ (g + 5(-))sgnK) (W^((7„) + -W"{Ur,] 



= J2 [2PW - (P + P(+)) {W{U:) + -W"iU+) 
where 2t/+ = u„ + u+. SinceEj 



LU„>0 



(B12) 



Y, [P-P{-)) [wiUn) + -W"{Un) )sgn{nn) 



W{U^) + -W"{U^) = -UnW'{Un), 



(B13) 
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the contribution of the last term to the final results is at least of order co^. The remaining terms give, after expansion 
in the external frequency, 



r^[ijsgiiK)H'(«„) - 5(9 + 9(-))sgn(n,j(w.'(CI„) + ^W'lUr., 



-^E [»'"("-) + (;:^)('*"'-' + 5»""<-' 



w„>0 

1 V sin 9 
27ri 2A 



I VZsin6'y Vi;sin6l/L KlsineJ 2 KlsmeJl) 



(B14) 



APPENDIX C: FLUCTUATIONS OF THE ORDER PARAMETER. 



Our starting point here is Eq. ( 102 ) for the linear correction to the anomalous propagator 



/i = 



evA(/ -/(-)) 



ILVo 



+ Kg + .9(-)) 2r2„ + v(V - 2ieA) 



Ai 



+i{2T)-\g + gi-)) 2r!„ + v(V - 2ie^) (A) 



-*(2r) 



-ievA(g-5(-)) 



lujo 



2n„+v{W~2ieA) (/) + (/(-)) 



If we define the angular average of /i by 

(/)=5(a-„)|l), 
we find, after carrying out the angular integration and ignoring terms of order A/l 



V2 



lUjQ 



(CI) 



(C2) 



5 = V^- r^htevAsinOl - A=) ^^' ^i ^^'^ ^ + i(<? + <?(-))2Csgn(aO (w{Un) + W'{U,,)) I. (C3) 



We now use Eq. ( |103| ) to determine the amplitude of the excited mode of the order parameter C = TTgN{0) J2n ^ 
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w „ Jo 2 V 2 
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/O 'i 1^0 

We can now use the gap equation to eliminate the need for a frequency cutoff and obtain 
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'V / d0\igsgniL,„)W{ 
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It follows that to leading order 
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We a lso g ive the expression for the distribution function /i. Neglecting the contributions of order {A/l), we use the 
Eq. (All) to compute 



28 



/i=evA -^ /^ ' +i(g + g{-))\2nn + v(V-2ieA)] \i (C7) 

= ^A/^eAylcos0 V ^e™*f ^"j™ x i ^W^^^^sgn^+iK) - 5(-)M^(")(-)sgn'"+iK-) llm) 

m— 

The evaluation of C is analogous to the calculation given above. 

In the transverse response calculation our starting point here is Eq. ( [123[ ) for the linear, in cyclotron frequency, 
correction to the anomalous propagator 



6f = [2f7„ + v/(V - 2ieA)] '|-^a(v - 2ieA) (/ - /(-)) + tASg, + iSA, (g + ff(-)) j- 



(C8) 



The solution of this equation follows exactly the steps described in the previous section. First we solve for the 
coefficient SC in SAi = SC\l). As 5C = J2n J '^^^^fi the denominator of the expression for SC is the propagator for 
the first excited mode of the order parameter, as it was for C. To evaluate contribution of each of the driving terms 
we notice that for our choice of A 

A^V - 2ieA) = AV^ = A^- a-a^ . (C9) 

Explicit evaluation of this term using the expression for the function / from Eq. (p^) shows that it contributes at 
order A/l compared to leading order terms. The driving term iASge only contributes in the intermediate frequency 
range since 

dge = -2iLUcevA sin sin (j)^--S-^ , (CIO) 

and, to the order in which we work, g — g{~) vanishes in the outer range. Then in analogy to the solution outlined 
above we obtain 

SCi = 6Ci = ieAA^f2[ujT)[uj^T). (Cll) 

and, for the angular average of the function 5/ needed to calculate the Thompson contribution to the conductivity, 

{bh) = V2^eA^^^ ^'^ f {sin 0(5 - g{-))W\U^) + 2i6c(g ^ g{-))^^x,{^r.){w{Ur.) + \w"{U,S)\\\) (C12) 

(<5/|> - y2^eA^^^^"^|sin0(.g-.g(-))vK'(t/„) + 2*,5c(5 + .g(-))sgn(17„)(M/(t/„) + \rw" {V,S)\^{\\. (C13) 

In the intermediate region g + 5(— ) = 0, and only the first term in each function contributes to the conductivity to 
leading order in (A//). 
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FIG. 1. Longitudinal resistivity as a function of the reduced magnetic field. 

FIG. 2. Transverse conductivity as a function of the reduced magnetic field. 

FIG. 3. Hall angle as a function of the reduced magnetic field. 
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